Full counting statistics of a charge pump in the Coulomb blockade regime 
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We study full charge counting statistics (FCCS) of a charge pump based on a nearly open single 
electron transistor. The problem is mapped onto an exactly soluble problem of a non-equilibrium 
g — 1/2 Luttinger liquid with an impurity. We obtain an analytic expression for the generating 
function of the transmitted charge for an arbitrary pumping strength. Although this model contains 
fractionally charged excitations only integer transmitted charges can be observed. In the weak 
pumping limit FCCS correspond to a Poissonian transmission of particles with charge e* = e/2 
from which all events with odd numbers of transferred particles are excluded. 



c3 



i 

-a 
c 

o 
o 



(N 
> 



O 



i 

C 

o 
o 



X 



PACS numbers: PACS numbers: 73.23.Hk, 73.40.Gk, 72.10.Bg, 71.10.Pm 



Charge pumping has attracted considerable theoretical 
and experimental interest. It occurs when the Hamilto- 
nian of the system changes slowly in comparison with 
the characteristic time scales of the problem. At the 
end of the pumping cycle, when the Hamiltonian re- 
turns to its initial value, a finite charge may be trans- 
mitted through the system. The amount of the trans- 
ferred charge depends on the details of the pumping cy- 
cle. This idea is due to Thouless,EJ who showed that in 
certain one-dimensional systems the transmitted charge 
is quantized in the adiabatic limit. Most of research ef- 
forts have focu^djj^chjs^ through meso- 
scopic devices.ra$™$BBffl 

Motivated by the efforts to build an accurate standard 
of electric current most experiments concentrated on sin- 
gle electron pumps in which the charge pumped during 
one G-yeJe is quantized due to the Coulomb blockade ef- 
fectsfltl Such devices are already used in metrological apt 
plications to produce an accurate capacitance standard.EI 

Understanding of noise properties of the pumped cur- 
rent and of the accuracy of quantization of the pumped 
charge are very important for metrological applications. 
In this case it is desirable to know not only the av- 
erage pumping current and its second moment (noise 
power) but the whole distribution function of the pumped 
charge. Such full charge counting statistics (FCCS) were 
first considered in Refs. [l4|,[l5| for systems with non- 
interacting electrons. 

In the present paper FCCS for a charge pump based 
on a single electron transistor are considered. More pre- 
cisely, the device in question consists of a quantum dot 
connected to the left and right leads by single channel 
quantum point contacts labeled by the index a = ±1, 
see Fig. [l]. Such devices, can be fabricated in semicon- 
ductor heterostructuresu where the electrons in the two- 
dimensional electron gas (2DEG) in a heterostrocture are 
electrostatically confined to the area of the dot by a nega- 
tive voltage which is applied to the metallic gates located 
on top of the 2DEG. The reflection amplitudes r a in the 
contacts are controlled by the voltages on gates a and 
are assumed to be small throughout the pumping cycle, 




FIG. 1: Schematic drawing of a single electron transistor 
electrostatically defined on a surface of a two dimensional 
electron gas. The quantum dot is connected to two leads by 
single channel quantum point contacts labeled by a. The 
voltages on the gates G and ±1 determine respectively the 
average electron number in the dot, N(t) and the reflection 
amplitudes, r±i(t) in QPCs. 



r a Cl. The Coulomb interaction of electrons in the dot 
can be treated within the constant interaction model: 



H c =E c (N-N(t)Y 



(1) 



where N is the number of electrons in the dot, Ec is the 
charging energy and N(t) is the dimensionless parameter 
proportional to the voltage on the central gate G. 

At low temperatures, T <C E Cl the electron transport 
across the device is dominated by cotunneling processes. 
The quantum dot is assumecL-tp be sufficiently large so 
that elastic cotunneling effectstll can be neglected and the 
transport of electrons across ±he device is dominated by 
the inelasting cotunneling]i§E3 In addition the electrons 
are assumed to be spin-polarized. This can be realized 
experimentally by applying a strong magnetic field par- 
allel to the plane of the 2DEG. 

Sufficiently strong pumping can lead to a non- 
equilibrium distribution of electrons in the dot. Below 
we assume that the deviations from equilibrium may be 
neglected. This requirement imposes a limitation on the 
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number of pumping cycles in the absence of energy re- 
laxation in the dot. Indeed, inelastic cotunneling can be 
thought of as a coherent process in which an electron, 
say from the left lead, enters a certain quantum state 
in the dot and an electron from a different state in the 
dot leaves into the right lead. As a result of such a pro- 
cess an electron is transferred across the device and an 
electron- hole pair is created in the dot. Upon completion 
of N pumping cycles roughly N electron-hole pairs will 
be created. The number of electron-hole pairs in equi- 
librium may be estimated as T/5±, where Si is the single 
particle mean level spacing in the dot. Therefore, for the 
deviations from equilibrium to be small we assume that 
the dot is sufficiently large so that the mean level spacing 
S\ <C T/N. In the presence of energy relaxation in the 
dot this condition may be relaxed. 

At frequencies below the charging energy Eq the 
puniDuig cycle is described by a single complex vari- 
ablejjy z(t) = n(t) exp{iTrN(t)) + r_i(t) exp(-inN(t)). 
The average pumping current for this cycle was obtained 
in Ref. | 

Here we study FCCS for this pump. The probability 
distribution function Pn(Q) for the charge Q transmitted 
through the dot upon completion of N pumping cycles is 
determined by the generating function 



Fv(A) = ^exp(iAQ)P A r(Q), 



(2) 



where the charge Q is measured in units of the abso- 
lute value of the electron charge, e, and the sum goes 
over all its possible values. The n-th cumulant ((Q n )) of 
the transmitted charge may be determined from Fpf(X) 
through the relation 



((Q n )) 



d n hiF N {\) 



i n d\ n 



(3) 
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Below we concentrate on the pumping cycle in which 
z(t) = zoexp(— iujt). In this case we obtain the following 
generating function 

^ oo 

\nF N {\) = - In (1 + cos 2 0(e ; )[n-(£;)[l -«+(£;)] 

1 — — OQ 

x (e lA - 1) + n+(ej)[l - n_(ei)] (e~ a - l)]) - iWA.(4) 

Here e; = lo(21+1)/(2N), with I being an integer, denotes 
the discrete fermionic frequency, n±(e) — no(e ± lo) = 
(e' e+a ')/ T + l)~ 1 is the Fermi distribution function shifted 
by ±oj, and exp[£0(e)] = (e + iT)/y/e 2 + T 2 , where T = 
2-/\z \ 2 E c /n 2 , 7 = expC, with In 7 = C w 0.5772... 
being the Euler constant. 

Note that since the generating function F/v(A) in 
Eq. is periodic in A with the period 2tt, Fn(X + 2tt) = 
Fzv(A), only integer values of charge Q can be transmit- 
ted. 

Eq. (|4|) acquires a particularly simple form at low tem- 
peratures, T <C lo. Approximating the Fermi functions 



by the step functions n±(ei) — 8(— ei =p uS) we observe 
from Eq. (^) that only the energy interval — lo < ei < lo 
contributes to the pumped charge. Using the Poisson 
summation formula we can write Eq. (^) as 



l/T 



lnFjv(A) = -iNX + TN 



n= — oo 



l/T 



xln.x 2 e lA/2 + l-lna; 2 + l, 



(5) 



where T = T/lj is the relative pumping strength. For 
long observation times tT — 2irYN ^> 1 the terms n^O 
become small due to the presence of quickly oscillating 
factors in their integrands. Explicit evaluation of the 
main, n — 0, term gives 



In F w (A) = AT In 



T 2 e" lA + 1 
T 2 + l 



2TiVarctan(T- 1 ) 



-2T7Ve~ lA/2 arctan ^T -1 e a/2 



(6) 



With the aidofEq. (|) we obtain for the average pumping 
current, 



I = 



£M = _fE arctan(T -i). 

r 2tt 



(7) 



The initial growth of the current with the pumping fre- 
quency lo saturates at / = — eT/4 for large lo. 

In the strong pumping limit, T>1, Eq. dq) yields 
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In F w (A) = -iNX+ ^ (e jA - 1 



(8) 



where the first term contributes only to the average cur- 
rent, and the second term describes a Poisson process for 
particles with an integer charge e* = e and transmission 
frequency 10 j (67rT 2 ). 

In the limit of weak pumping, T <C 1, we can per- 
form the integration over x in Eq. (JsJ) from — 00 to 00. 
Retaining the terms with n/0we obtain 



F w (A) 



cosh [7rT7Ve- iA / 2 ] 
cosh [ttTTV] 



(9) 



To evaluate the cumulants (||) it suffices to know Fjy(A) 
at A — > 0. For long observation times rr > 1 we can 
neglect the exponentially small terms in Eq. (|^) and write 
the logarithm of the generating function as, 



In Fat (A) = irYN (V iA / 2 - l) 



(10) 



This formula can also be derived directly from Eq. (g). 
It corresponds to a Poisson process which describes in- 
dependent transmission of quasiparticles with the av- 
erage transmission frequency T/2 and fractional charge 
e* = e/2. The true limiting expression for weak pump- 
ing, Eq. @ is periodic in A with the period 27r, allow- 
ing transmission of only integer charges. It is easy to 
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check that Eq. describes a Poisson process for charge 
e /2 particles from which all transmission events with odd 
numbers of transferred particles have been excluded. The 
corrections to Eq. (||) are small ~ T 2 . 

One may define the effective charge e* of the carri- 
ers through the ratio of the variance of the transmit- 
ted charge to its average value for intermediate pump- 
ing strengths as well. However Eqs. show it can 
be interpreted as a charge of independently transmitted 
particles only in the limits of weak, T < 1, and strong, 
T > 1, pumping. The coefficients in the Taylor expan- 
sion of Eq. (^) in powers 1/T™ (for strong pumping) or 
in powers T n (for weak pumping) represent Poissonian 
transmission processes of multiple charge ne*, in agree- 
ment with Ref. 

Below we present the derivation of the above results. 
At T,lu <C Ef£H the pumping cycle is described by the 
Hamiltoniarj 



H 



dk 



2tt 



-z*(t) 
z(t) 



ill) 



where k = yj ^vEq j^" 1 ■ In Eq. (|Tl]) tyf. is a vector fermion 
operator in Gorkov-Nambu notations and is expressed 
through the creation and annihilation operators c k and 
c\ as fy k = (c^,c-fc), and 173 is the Pauli matrix. In 
this model electrons have a linear spectrum e& = vk and 
are coupled to a resonant state described by a Majorana 
fermion £, Q 2 = 1. The current through the pump is 
given by 



I = - 



evp 



dk 



(12) 



For the pumping cycle considered here the gauge trans- 
formation 



*fe — > exp(icr 3 wi)*fc, 



(13) 



removes the time dependence of the Hamiltonian, z{t) — > 
Zq. As a result, the chemical potentials of electrons and 
holes shift by ±w. The current operator in this gauge ac- 
quires an additional anomalous term and takes the form, 



1 = -— dk*l* k + -. 



(14) 



The stationary Hamiltonian ( |il| ) (with z(t) — > z ) was 
diagonalized by MatveevEJ in terms of the linear combi- 
nations of particle and hole operators, 



Ci. — 



v/2 ' 
e k ± iT c k 



(15a) 



VfT 



r f de k > 



Ck' 



-(15b) 



Both signs in Eq. (15b) give equivalent expre ssion s. For 
the upper/lower sign the last term in Eq. ( |l5b| ) gives 
a vanishing contribution to ^(x) at x — > ±00 after a 
Fourier tran sform ation to the real space. (The second 
term in Eq. (15b) corresponds to the resonant state van- 
ishing for x — > ±00.) 

Having observed these asymptotic properties of oper- 
ators ( |l5|) we can readily build scattering states corre- 
sponding to the scattering of an electron, 




e- l9 ^Ck, 



-00, 



Ck cos8(ek) - ic_ k smd(ek) , x — > +00, 



(16) 



and a hole 



e -ie(e k ) Ck _ Ck 

p -i9{e k ) t 



X —> —OO, 



— ick sin0(efc) + c\ cos8(e k ) , x — * +00 



'(17) 



where cos8(e k ) = e^/y'^+T 2 , and sin0(efc) = 
V € 'k + T 2 - We can now write the scattering matrix 



cosO(ek) —ismO(ek) 
-ism9(e k ) cos9(ek) 



(18) 



for the scattering between electron (positive current) and 
hole (negative current) states. 

Thus, the problem reduces to the problem of non- 
equilibrium chiral fermions scattering off a resonant state 
at zero energy. The electrons and holes here are charac- 
terized by non-equilibrium distribution functions n±(e) 
defined below Eq. (|J) and may be represented by the 
diagonal matrix 



n(e) 







n+(e) 



(19) 



The full counting statistics for non-equilibriurp_non- 
interacting fermions have been extensively studiedo The 
generating function of the transmitted charge is given by 

F N (X) = exp(-iAA) x 

exp | Tr £ ln I 1 + ™( e *) (sl x (tk)Sx(e k ) - l)] jap) 

where S±\(e k ) = exp(±|<7 3 A)S(efc) exp(±|cr 3 A) with 
S(ek) defined in Eq. (|l8|). The first term in this equation 
arises from the anomalous term in the current operator, 
Eq. (111). Since electron and hole operators describe the 
same physical states, the sum over energies is restricted 
to positive frequencies in order to avoid double coun ting 
of degrees of freedom. Substituting Eqs. (|l^) and ( |T9| ) 
into Eq. rtgfj) we obtain the generating function (0). 
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In conclusion, we have obtained full counting statistics 
for a charge pump based on a nearly open single elec- 
tron transistor. In the spin-polarized case the problem 
is mapped onto an exactly soluble chiral fermion model, 
Eq. (11). In the weak pumping regime T <C u> the gen- 
crating function given by Eq. (|9|) corresponds to a Pois- 
son process of charge e* = e/2 particles with transmis- 
sion rate T/2 from which all events with odd numbers of 
transferred particles are excluded. Alhough all the mo- 
ments of the transferred charge obtained from Eq. (||) 
are practically indistinguishable from those of a simple 
Poisson process for charge e* — e/2 particles only inte- 
ger transferred charges may be observed in a pumping 
experiment. Since the Hamiltonian ( fil| ) of this model 
describes aj=| Luttinger liquid with an impurity one 
may expect that similar conclusion hold for other cou- 
pling strengths g ^ \ and other problems with fraction- 
ally charged excitaiipiis_realized for example in Quantum 
Hall experiments.BBEl 

Equation differs from the weak pumping result in 
Rcf. 13. The reason for this discrepancy and for charge 
fractionalization lies in the failure of perturbation the- 
ory in the reflection amplitudes r±± for our model at 
sufficiently low energies e < T. The effective reflection 
coefficient T 2 /(e 2 + T 2 ) that determines the strength of 
pumping approaches unity in this energy range. Thus, 
the true expansion parameter at weak pumping is not 
the reflection amplitude r±\ but the ratio of energy scales 



T = T/ui. Perturbation theory in the reflection ampli- 
tude fails for Luttinger models with an impurity at other 
interaction strengths g ^ \ as well wich leads to the ap- 
pearance of an energy scale analogous to Ilbelow which 
the system is in the strong coupling limit .E3 Ell 

Although we have focused on the low temperature case 
r«w the validity of the result (^) is restricted only by 
the condition T < Ec- Using Eq. (|J) we obtain the 
general expression for the average pumping current, 



eto 
2^ 



e 

47T 



e 2 de 



sinh^ 



e 2 +T 2 cosh^ 



cosh ■ 



(21) 



This formula reduces to the result of Ref. || in the linear 
response regime u> <C T. 

The case of zero pumping and finite external bias V 
can be obtained from the above expressions by substi- 
tuting lj — ► eV and omitting the anomalous term in the 
current. For example the non- linear I — V characteristic 
is obtained in this way from the second term of Eq. (|2l]) . 
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